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I give an outline of my recent proposal to take the QCD functional determinant in lattice simulations partially into account: The 
determinant is split into two factors, the factor referring to a standard background in each topological sector is kept exactly, the 
factor describing the effect of the smooth deviation of the actual configuration from the reference background is replaced by one. 
The issue of how to choose the reference configurations is discussed and it is argued that "topologically unquenched QCD" is an 
interesting starting point to study full QCD in lattice simulations as it gets the main qualitative features right from the beginning. 



1 What's wrong with quenched QCD ? 



2 What is "topologically unquenched QCD" ? 
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In QCD the fermion functional determinant is a nonlocal 
contribution to the gluon effective action. This nonlo- 
cality has an unpleasant effect in lattice calculations as 
it slows down present numerical algorithms dramatically 
when the quark-masses get small. 

In order to elude this problem most numerical simu- 
lations in the past have been done in the "quenched ap- 
proximation" where the determinant is simply replaced 
by one Id. More recently, simulations in the "partially 
quenched approximation" have become availableH, where 
the determinant gets evaluated at quark masses which 
are higher than those in the propagators. Thus (partial) 
quenching amounts to suppressing the contribution of all 
internal fermion loops in QCD by giving the quarks un- 
phenomenologically high or infinite masses. 

Attempts to introduce the corresponding modifica- 
tions in the low-energy theory artificially in order to learn 
how to correct for them - the results being "quenched" 
and "partially quenched Chiral Perturbation Theory" - 
have shown that the (partially) quenched approximation 
is, in some aspects, fundamentally different from the full 
theory: First of all, numerical results won in quenched or 
partially quenched simulations should (at least in princi- 
ple) be corrected Jbr the occurrence of "enhanced chi- 
ral logarithms" trU In addition, the if was found to 
be a pseudo-Goldstone boson in quenched QCD (as op- 
posed to the situation in QCD, where it is heavier than 
the lowest-lying octet of pseudo-scalar pseudo-Goldstone 
mesons by more than a factor y3) and its propagator 
shows - in case the low-energy analysis is correct - a 
pole of order two (which spoils any field-theoretic in- 
terpretation) right at the same position in the p 2 -plane 
where its first-order pole is Ho. Even without particle- 
interpretation, quenched QCD is strictly confined to Eu- 
clidean space-time; there is no continuation of its Green's 
functions to Minkowski space-timela. 

Thus it seems legitimate to search for an alternative 
starting point for approaching full QCD which gets some 
qualitative aspects of full QCD right from the beginning. 



The aim is to identify a part of the functional determi- 
nant which is cheap from the computational point of view 
but essential from the field-theoretical point of view. 

We start from the generating functional of (eucli- 
dean) QCD in the form where the fermionic degrees of 
freedom have been integrated out 
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det{fl + M) 
det(^+M) 



-GG 



(1) 



igA^ 
1 . 



2 ^ixuapG ' up 



is the (euclidean) Dirac op- 
the dual of the field- 
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erator and 

strength operator and where the measure DA includes 
gauge-fixing and Faddeev Popov terms. In ([j]) a factor 
which does not depend on the gauge field to be inte- 
grated over has been pulled out of the normalizing fac- 
tor and the convention is that the quark mass matrix 
M is diagonal and of rank Nf with the CP-violation 
stemming entirely from 9 (if 6 ^ ttZ); the shorthand- 
notation to be used here and in subsequent formulas is 

fNf 



Y[i=i det(p + irii) in the determinant and 
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detQp + M) 

^{IJH-My 1 ^ = YsiJiV^iP+rni) -1 ^) m the propagator. 

QCD is known to show a topological structured: The 
S'?7(3)-gauge-field configurations on R 4 with finite ac- 
tion boundary condition or on the torus T 4 fall into in- 
equivalent topological classes (labeled by an index v — 
g 2 /32iT 2 ■ J GGdx £ Z). In a given sector v, any two 
configurations may be continuously deformed into each 
other but not into any configuration with a different in- 
dex v. Due to this topological structure, the integral 
in (|l|) may be rewritten as a sum of integrals over the 
individual sectors J DA — > X^ez / DA^ and the deter- 
minant factorizes (If) = p^') 
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det(^+M) det(#+A/) 
where the first factor depends on v only and thus may 
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be pulled out of the integral. In this form it is obvious 
that quenched QCD actually does two modifications: It 
sets both determinant factors equal one, whereas par- 
tially quenched QCD keeps both of them at the price of 
using unrealistically high quark masses. 

In "topologically unquenched QCD" the first factor 
in (J2j> is kept exact (with the same quark mass as in the 
propagator between rj, n) and only the second factor is 
replaced by oneQ, i.e. the "theory" is defined through 

J DAM eViP^+My'v e -Ji GG . ( 3 ) 

At this point, the motivation to treat the two deter- 
minant factors in (|^) on unequal footing is just an eco- 
nomic one: The "topological" factor is universal for all 
configurations within one class and bears the knowledge 
about the nontrivial topological structure of QCD. On 
the other hand, the "continuous" factor in (||) causes a 
dramatic slowdown in numerical simulations as this part 
of the determinant (or its change) has to be computed 
for each configuration individually. 

Note that, in contrast to the full-QCD generat- 
ing functional (|l|), the "topologically unquenched" trun- 
cation (H) does depend on the choice of reference- 
backgrounds. Later, two alternative strategies of how to 
choose, in a given sector, the reference-configuration are 
described and such a choice of strategy, once it is done, is 
considered part of the definition of the theory. Obviously 
there is no a-priori evidence that - with any of these two 
choices - the approximation (^|) should be particularly 
good. Nevertheless, known analytical knowledge about 
QCD in a finite box seems to indicate that including the 
"topological" part of the determinant is sufficient to get 
some basic features of full QCD qualitatively right. 

3 How to implement "topologically unquenched 
QCD" ? 

Obviously, for "topologically unquenched QCD" to be 
practically useful, at least two requirements have to be 
fulfilled: First, a recipe of how to make a good choice for 
the reference backgrounds in (||) is needed. Second, the 
costs in terms of CPU-time for separating the "topolog- 
ical" part of the determinant from the remainder have 
to be smaller than the costs would be to compute the 
determinant as a whole. These issues shall be discussed. 

3.1 Upgrading a quenched sample 

A quenched sample may get modified to be representative 
in the sense of "topologically unquenched QCD" : 



1. Use a method you consider both trustworthy and 
efficient to compute for each configuration its topo- 
logical index v. 

2. Use the gauge action you trust to compute, in 
each class, the gauge-action of every configuration 
as well as the class-average and choose the 
reference-configuration according to one of the fol- 
lowing two prescriptions: 

(i) Choose - out of the class v - the configuration 
with minimal gauge-action as the representa- 

(ii) Choose - out of the class v - the configuration 
for which the gauge-action is closest to the 
class-average as the representative A^fy. 

3. Use the fermion action and the method you con- 
sider both trustworthy and efficient to compute the 
determinants det((^]+M)/((^> +M)). 

4. For any higher topological sector either include the 
corresponding determinant computed in step 3 into 
the measurement or eliminate the corresponding 
fraction of configurations from that sector. 

A few points need immediate clarification: 

The first problem is that on the lattice, the topo- 
logical structure of the continuum-theory is washed out; 
simply computing v = g 2 /32tt 2 ■ J GG dx gives a value for 
v which is, in general, not an integer. A solution to this 
problem is so crucial to the overall performance of the 
"topologically unquenched" approximation that it shall 
be discussed in a separate subsection. 

The second point is that for computing the determi- 
nant ratio in step 3 one cannot rely on any method which 
is tantamount to an expansion in 5 A, since the two back- 
grounds are far from each other. The necessary ab-initio 
computation is achieved e.g. by the eigenvalue method: 
typically the first few hundred eigenvalues of the Dirac 
operator on a given background may be determined. 

Finally: What's the difference in terms of physics be- 
tween the two strategies of how to choose the reference- 
backgrounds in step 2 ? We emphasize that for either 
choice there is a sound theoretical motivation. Strat- 
egy (*) ~~ choose the configuration which minimizes the 
gauge action - is nothing but the semiclassical ansatz 
being pushed to account for topology: Within each sec- 
tor, the determinant is exact for the configuration having 
least gauge-action, i.e. for the one which, in a semiclas- 
sical treatment, gives the dominant contribution of that 
sector to the path-integral. Strategy (ii) - choose the 
configuration which, in its gauge-action, is closest to the 
class-average of that sector - takes into account that the 
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Monte Carlo simulation as a whole doesn't try to mini- 
mize the total action density but rather the free-energy 
density: The configuration which is most typical in a cer- 
tain sector is not the one with minimal action but the one 
which has additional instanton-antiinstanton pairs plus 
topologically trivial excitations such as to find an opti- 
mum between the additional amount of action to be paid 
and the additional amount of entropy to be gained. It is 
the very aim of the second strategy to choose in each sec- 
tor a "most-typical" background (which realizes such an 
optimum pay-off) as reference-configuration. It should 
be stressed that even though the two strategies end up 
selecting reference-backgrounds which look highly differ- 
ent (exceedingly smooth in the first case versus pretty 
rough in the second case) the final results may still be 
close to each other - the only thing which matters is the 
(strategy-intrinsic) determinant ratio computed in step 3. 

3.2 Generating a "topologically unquenched" sample 

The fact that the "topological" determinant in @ is 
a number which depends only on the total topologi- 
cal charge of the actual configuration but not on its 
other details suggests that one could try to precom- 
pute these "topological" determinants on artificially con- 
structed backgrounds prior to running the simulation. 

Within the strictly semiclassical strategy which is re- 
flected by choice (i) the reference-backgrounds are got- 
ten in a rather simple way: Place v copies (for v > 0) 
of a single-instanton solution with typical radius (i.e. 
p ~ 0.3fm, cf. H) on the lattice and minimize the ac- 
tion with respect to variation of their relative orienta- 
tions and positions. Within the more realistic strategy 
of choice (ii) which accounts for the competing effects of 
increased action versus increased entropy the reference- 
backgrounds are constructed as follows: Place v instan- 
tons (for v > 0) with typical sizes (i.e. p ~ 0.3fm, 
cf. a) randomly on the lattice plus additional instanton- 
antiinstanton pairs such as to achieve a total instanton 
density of lfm~ 4 (cf.B). Optionally, this background may 
be dressed with thermal fluctuations by applying a rea- 
sonable number of heating-steps (monitoring v in order 
to guarantee that it stays unchanged). 

Having defined the standard backgrounds in this way 
a pure Metropolis algorithm generating a "topologically 
unquenched" sample employs the following steps: 

1. Make use of an updating procedure to propose a 
new configuration and determine its topological in- 
dex (via the method you trust). 

2. If the configuration realizes a previously unseen v: 
Evaluate the functional determinant on the stan- 
dard background constructed for that v according 
to your choice of strategy. 



3. Base the decision on whether to accept the pro- 
posed configuration on 
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denotes the gluonic action of the 



newly proposed / last accepted configuration. 

There is one conceptual deficit this algorithm suffers 
from: The procedure for constructing the standard- 
backgrounds makes use of knowledge about the size- 
distribution and partly about the density of (anti-) in- 
stantons which was won in previous lattice-studies. In 
other words: The "topologically unquenched" simula- 
tion as outlined above is not entirely from first princi- 
ples. Moreover, the quality of the final sample depends 
on how appropriate the artificial backgrounds are which 
were used in the computation of the standard deter- 
minants. Constructing these reference backgrounds is 
particularly demanding within strategy (ii) as it means 
that one has to do an a-priori guess which configuration, 
within a given sector, is "most-typical" in the sense of full 
QCD. Even if "most typical" is translated into a technical 
criterion (e.g. the configuration which, in its total effec- 
tive action is closest to the corresponding class-average 
of a finite full-QCD sample) there is no other algorithmic 
solution to this problem than by doing a full-QCD sim- 
ulation. Thus it is reasonable to construct the reference 
backgrounds as indicated above, thereby making use of 
existing knowledge concerning size-distribution and den- 
sity of instantons in QCD - in particular as there is, for 
strategy (ii), a final check concerning the quality of the 
artificial backgrounds: In case the guess would have been 
just perfect, the a-priori ratios of determinants evaluated 
on these backgrounds (i.e. the ratios which got used in 
the "topologically unquenched" simulation) would per- 
fectly agree with the analogous a-posteriori determinant 
ratios evaluated on the "most-typical" configurations as 
produced by the run. Accordingly, a procedure one might 
think of trying in case the agreement turns out to be less 
than satisfactory is just to start over with the simulation 
- but this time using the "most-typical" configurations 
found in the first run rather than the artificial guesses. 

3.3 Measuring topological indies 

Determining the topological index of a newly proposed 
background in a way which is fast and reliable is so crucial 
to the overall-performance of "topologically unquenched 
QCD" as to justify few remarks about this point. 

There are several methodsEI to determine, for a given 
configuration, its topological index v. Some of them were 
recently compared and found to give - when implemented 
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with sufficient care -_to comparable results for the topo- 
logical susceptibility E9. Nevertheless, it is clear that in 
the present case, where nothing is known about the spec- 
trum of the Dirac operator on the background at hand, 
the field-theoretic method is likely to determine v in the 
quickest possible way. However, the fact that on the 
lattice the relevant operator undergoes thermal renor- 
malization provides a challenge: Simply integrating the 
Chern density, i.e. computing g 2 /(327r 2 ) f G a G a dx 
gives a value which is, in general, not close to an inte- 
ger; in fact, a histogram-plot over many configurations 
tends to reveal accumulations near regularly displaced, 
non-integer values, e.g. near 0, ±0.7, ±1.4 etc. There are 
two options of how to deal with this situation: 

The first, simplistic, approach is just to define a "con- 
fidence interval" - e.g. ±0.2 - around each of the values 
0, ±0.7, ±1.4, . . . and to assign the configurations lying 
within these bounds the indices v — 0, ±1,±2... etc. 
The remaining configurations which didn't get an index 
assigned are then simply tossed away. 

The second, more sophisticated, approach is to make 
use of the fact that cooling a configuration is able to re- 
move the effect brought in by thermal renormalization: 
Cooling a set of gluon-configurations results in the peaks 
(in the histogram plot) being shifted closer to the cor- 
responding integers and the valleys between the peaks 
getting thinned out under each sweep. 

The problem is that these two methods do not nec- 
essarily agree in their results for a given configuration - a 
fact which can be understood on rather simple grounds: 
Under repeated cooling with the naive (Wilson) action, 
a single-instanton solution shrinks monotonically until it 
finally falls through the grid. In order to prevent the cool- 
ing algorithm at least from loosing the large instantons 
one has to modify the action w.r.t. which cooling is done 
in such a way that all instantons with a radius p above 
a certain pthr tend to get blown up ("over-improved ac- 
tion") or stay constant ("perfect action") under a sweep, 
where typically ptnr — 2.3a. From the evidence given 
in EH how quickly cooling with an "over-improved" ac- 
tion tends to pin down g 2 / '(327r 2 ) f G^G^ dx near 
an integer (say 5 sweeps to be within 2.99 and 3.01, 
etc.) performing 0(3) "over-improved" cooling sweeps 
seems to be sufficient to get an unambiguous assignment. 
The price to pay, however, is that the small instantons 
(p < Pthr) get compressed and finally pushed through 
the grid even more efficiently than under cooling with an 
unimproved action til. From these consideration we con- 
clude that the field-theoretic method with cooling yields, 
once it has stabilized, a correct assignment for the lat- 
ter cooled configuration which, however, isn't necessarily 
appropriate for the initial configuration which may have 
contained small (p < 2.3a) instantons. On the other 
hand, determining the topological index by the first (sim- 



plistic) approach (no cooling being involved) has an infe- 
rior signal-to-noise ratio (about half of the configurations 
can't get assigned an index and have to be tossed away) 
but for the remaining ones the procedure is sensitive to 
all instantons the lattice can support (i.e. p > 0.7a). 

3.4 Rudimentary cost analysis 

The overhead as compared to a quenched simulation re- 
sults from the CPU-time spent on determining v for ev- 
ery newly proposed configuration and from the deter- 
minants which get evaluated. Preparing the reference- 
backgrounds and computing the determinants is a fixed 
investment which is given by L,a,m only (i.e. indepen- 
dent of the length of the simulation) and evaluating the 
O(10) standard determinants (for_ nowadays typical val- 
ues of m and L) is pretty cheap 0. On the other hand, 
determining for each configuration its index v gives rise 
to costs which grow linearly in simulation-time and thus 
provide the main overhead (as compared to Q-QCD) in 
a long run. As a consequence, the method for determin- 
ing the topological index will have the greatest impact 
on the overall-performance in TU-QCD. We have advo- 
cated choosing a field-theoretical definition with little or 
no cooling at all, which means that either 0(3) cooling- 
sweeps are performed or 50% of the configurations have 
to be tossed away. Accordingly, in an approximation 
where a cooling-sweep is considered twice as expensive as 
a complete update, the overhead from ^-determinations 
is roughly a factor 2... 6 over a quenched simulation. Thus 
doing a "topologically unquenched" run might be consid- 
ered an alternative to a high-statistics quenched run. 

4 What about qualitative features of "topologi- 
cally unquenched QCD" ? 

In "topologically unquenched QCD" a determinant is in- 
troduced which - as is seen from (|j) - only influences the 
relative weight of the different topological sectors; within 
each sector there is no difference to quenched QCD. 

For QCD in a finite box Leutwyler and Smilga have 
shown that in the regime [] USm ^> 1 the distribution of 
topological indices is gaussian with width£3 

{v 2 ) = VEm/Nf . (5) 

In quenched QCD the corresponding distribution is much 
broader as there is no determinant which suppresses the 
higher sectors. In "topologically unquenched QCD" the 
standard determinants result in the higher sectors being 
suppressed as compared to a quenched sample but the 

a S = lim m ^o limv^oo K 1 /'' 1 /')!; where irij = ra Vi for simplicity; 
note that VEm — > 00 when m — > as the box has to be scaled 
accordingly: L ~ 1/M W ,M^ ~ ^hnd m - 
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amount of suppression strongly depends on the strategy 
for selecting or constructing the reference backgrounds. 

In strategy (i) a sectorial determinant is introduced 
which is exact for the background which - from the clas- 
sical point of view - dominates that sector. The point 
is that this semiclassical treatment is indeed justified for 
sufficiently small coupling-constant, i.e. in a ridiculously 
small box where the topological distribution in QCD is 
known to be extremely narrowE 2 ! As the box-volume in- 
creases the effective coupling gets stronger and strategy 
(i) is unable to account for this change. To see this more 
clearly we stipulate the validity of the index theorem on 
the lattice □ which allows us to rewrite the_two factors 
in pi) using the Vafa-Witten representationllj 
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Strategy (i) retains a determinant which is appropriate 
in a small volume and thus strongly suppresses the higher 
topological sectors. As it comes to larger volumes, the 
semiclassical treatment breaks down and the quantum 
fluctuations packed into the "continuous" determinant 
(0) prove able to milder the suppression - in full QCD, 
but not within strategy (i). The virtue of strategy (ii) 
is that this change is accounted for by successively re- 
defining the standard-backgrounds used in (^). In other 
words: Within strategy (ii) parts which would belong to 
(0) in (i) are gradually reshuffled into the "topological" 
part (H) as the box-volume increases. As a consequence, 
either strategy is supposed to be trustworthy as long as 
VT,m < 1, but only strategy (ii) may give a reasonable 
approximation to full QCD in the regime VT,m 1. 

Comparing the two factors (|^) and ([?]) one ends up 
realizing that the "topological" determinant (^) has ex- 
actly the same structure as its QCD counterpart (the lat- 
ter comes without the subscript "std" in the numerator) : 
The essential ingredient- is the prefactor m) v K In QCD, 
this prefactor is knownE3 to cause the strong suppression 
of nonzero indices in the limit VSm <§; 1 . The fact that it 
is still around in the "topologically unquenched" approxi- 
mation (with either choice for the reference-backgrounds) 
means that TU-QCD (unlike Q-QCD) shows the phe- 
nomenon of chiral symmetry restoration if the chiral limit 
is performed in a finite box. 

Finally, the fact that the number of virtual quark- 
loops is not restricted in TU-QCD means that there is 
an infinite number of diagrams contributing to the rf- 
propagator (not just the connected and the hairpin di- 
agram as in Q-QCD) and this propagator may even be 
well-defined in the field-theoretic sense. 



In summary, the fermions in "topologically un- 
quenched QCD" are fully dynamical, but they interact 
in a way which does not pay attention to the details of 
the gluon background configuration but to its topological 
index only and this seems to be sufficient to get a number 
of basic features of full QCD qualitatively right. 
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